(Appeared in Journal of Mathematical Physics 35, 1172-1376, 1994) 
Deformations of Gabor Frames 

Gerald Kaiser 

Department of Mathematical Sciences 
University of Massachusetts at Lowell 
Lowell, MA 01854, USA 
email: kaiserg@woods.ulowell.edu 

June, 1993 
ABSTRACT 

The quantum mechanical harmonic oscillator Hamiltonian H = (t 2 — df)/2 
generates a one-parameter unitary group W{9) = e l9H in L 2 (R) which ro- 
tates the time-frequency plane. In particular, W(tt/2) is the Fourier trans- 
form. When W(9) is applied to any frame of Gabor wavelets, the result 
is another such frame with identical frame bounds. Thus each Gabor frame 
gives rise to a one-parameter family of frames, which we call a deformation of 
the original. For example, beginning with the usual tight frame T of Gabor 
wavelets generated by a compactly supported window g(t) and parameter- 
ized by a regular lattice in the time-frequency plane, one obtains a family 
{Tq : < 9 < 2%} of frames generated by the non-compactly supported 
windows g e = W(9)g, parameterized by rotated versions of the original lat- 
tice. This gives a method for constructing tight frames of Gabor wavelets for 
which neither the window nor its Fourier transform have compact support. 
When 9 = 7r/2, Tq is the well-known Gabor frame generated by a window 
with compactly supported Fourier transform. The family \Tq\ therefore 
interpolates these two familiar examples. 

PACS numbers: 02.20.+b, 03.65.-w. 



1. Introduction 

If f(t) is a complex- valued, differentiable function, let 



(Qf)(t)=tf(t), (Pf)(t) = -if(t). (1) 

Q and P extend to unbounded, self-adjoint operators on L 2 (R) which satisfy 
the (Heisenberg) commutation relation [Q, P] = QP — PQ = i/, where J is 
the identity operator. Being self-adjoint, Q and P generate one-parameter 
unitary groups operating on L 2 (R), 

Ufa) = e iujQ , V(s) = e~ lsP . (2) 
U(u) and V(s) act by modulation and translation, respectively: 

{U{<")m) = f{t), (V(s)f)(t) = f(t - s). (3) 

The local (Lie-algebraic) commutation relation [Q, P] = il has the global 
(Lie-group) equivalent 

U(u)V(s) = e luJS V(s)U(u). (4) 
Hence the set of operators G(4>, u, s) = e 1 ^ U(uj) V(s) satisfies the relation 

G((f>, iv, s) G(cj)', J, s') = G((j) + (/)'- u's, u + u',s + s'), (5) 

so it forms a group of unitary operators on L 2 (R), known as (a representation 
of) the Weyl-Heisenberg group. The Gabor transform (or windowed Fourier 
transform) may be viewed entirely in terms of this group. (See Daubechies 
[1] or Kaiser [2] for general background on the windowed Fourier transform.) 
Namely, given a window function g(t) in £ 2 (R) with \\g\\ = 1, define 



g u , s = U(u) V(s)g. 
2 



(6) 



Then (3) shows that gu, iS (t) = e tU}t g(t — s), which gives the translated and 
modulated windows that form the basis of the continuous Gabor transform 
and its inverse: 

/oo 
dte-^g(t-s)f(t), 

(7) 



/(£) = -!- J j^duds guj,s{t) f{u,s). 



Under certain conditions, a discrete subset of such "Gabor wavelets" 
g u , s is sufficient to reconstruct f(t) from f{u, s). For example, suppose that 
g has compact support in [0, r] and choose an interval < T < r. Suppose 
that the periodic function 

H(t)^rJ2\9(t-nT)\ 2 (8) 
is bounded above and below by positive constants, i.e., 

< A < H{t) < B. (9) 
Then it can be shown (cf. Daubechies et al. [3], Kaiser [4]) that 

f(t)=H(t)- 1 J2 Y.9^s n {t)~f{^s n ), (10) 
where u m = 27rm/r and s n = nT. In fact, (9) is equivalent to 

^ii/ii 2 <EEk^-/)i 2 ^ s ii/ii 2 ' ( n ) 

which is just the condition that the set of vectors {g UJm ,s n '■ m,n G Z} form 
a frame in £ 2 (R) with frame bounds A, B. This frame can be made tight if 
we replace g by the new window 

h^^Hitr^git). (12) 
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Since h is also compactly supported, and T^ neZ \h(t — nT)\ 2 = 1, the above 
shows that {h UmySn } form a tight frame with bounds A' = B' = 1, i.e., a 
resolution of unity. Note, however, that \\h\\ 2 = T/t < 1, so this frame is 
an orthonormal basis if only if T = r. But if T = r, then g (and therefore 
also h) must be discontinuous in order to satisfy (9); this means that g(u>) 
and h(u) have slow decay, giving f(u, s) poor frequency resolution. (This is 
a special case of the Balian-Low theorem; see Daubechies [1].) 

The above construction depends crucially on the fact that g has compact 
support, since it expands g(t — s) f(t) in a Fourier series. A similar construc- 
tion can be made starting with a window which is compactly supported in 
frequency (expand in a Fourier series in the frequency domain). This leads 
again to a resolution of unity in terms of Gabor wavelets. For windows that 
are neither compactly supported in time nor in frequency, both constructions 
fail. For such windows, it is generally difficult to give explicit constructions 
of resolutions of unity. The theorem proved in the next section provides 
a method for generating a one-parameter family of Gabor frames starting 
from an arbitrary Gabor frame. If the original frame is tight or orthonormal, 
then so is each member of the family it generates. A general member of this 
family will have a window which is neither compactly supported in time nor 
in frequency, even when the original window has compact support. 

2. Deformations of Frames 

Consider the harmonic oscillator Hamiltonian 

H= 1 -(Q 2 + P 2 )= 1 -(t 2 ~d?), (13) 
which is an unbounded and self-adjoint operator on L 2 (R). Note that 

[H,Q] = -iP, [H,P] = iQ. (14) 
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Let W{6) = e l0H be the one-parameter group of unitary operators generated 
by H, and define the one-parameter families of operators 



Q(9) = W{9)QW{-9), P{6) = W{6)PW{-6). (15) 

Then 

±Q(0) = i [H, Q{9)] = i W{9) [H, Q] W{-9) = P(B), 

? (16) 
-P{9) = i [H, P(9)} = i W{9) [H, P] W{-9) = -Q{9). 

These operator equations (with the initial conditions Q(0) = Q, P(0) = P) 
can be integrated to give 

Q{9) = Qcos9 + Psm9, P{9) = Pcos9-Qsm9. (17) 

Hence, W(6) "rotates" the operators Q and P. W{9) acts on the "global" 
operators U{uj) and V(s) by 

W{9) U(u) = W{9) e lujQ = e lujQ{e) W(6), 

W{9) V(s) = W{6) e~ lsP = e ~ lsP ^ W{9). ^ 

From (17) and (18), we can compute the action of W(9) on the Gabor 
wavelets g^^: 

W{9) g Uj8 = W{9) e lujQ e~ lsP g = e tujQ(e) e" lsP(e) W{9) g. (19) 

To obtain an explicit expression, let g e = W{9) g. This is a new window 
function which generate a new family of Gabor wavelets g e w a = U{uj) V(s) g e . 

Theorem 1. W{9) acts on the Gabor wavelets g u , s by replacing the window 
g with g e ', rotating the labels (u>, s), and multiplying by a phase factor 7: 

W(9)g OJ , s = 1 (u,s,9)gi (e)M0) , (20) 



where 



and 



7(0;, s, 9) = exp 



- (u 2 - s 2 ) sin(20) + ius sin 2 9 



(21) 
(22) 



s(6>) = s cos 6> + cj sin 6>, cu(6>) = cj cos 9 — s sin 9. 
The new window g 9 is given by 

/CO 
dt exp [iut esc 9 - i(u 2 + t 2 ) cot 9] g(t) (23) 
-00 

for 9 ^ nir. When 9 = n/2, (23) reduces to g^ /2 {u) = e l7T ^g(u), where 
g is the Fourier transform of g. The singular cases are given by g n7T (t) = 
i n g((-l) n t), neZ. 



Proof: We need the identity 



e i(aQ+bP) _ e iab/2 e iaQ ^bP 



a, b G R . 



(24) 



which is closely related to (4). Setting a = cos 9 and (3 = sin6>, we have 



jQ(0) 



e iujaQ+iu>pP = exp ^2 a ^ 2 j e iu>aQ e iujf3P ^ 
e -isP(9) = jsfSQ-isaP = exp [_ zs 2 a/3/2 ] g - lsa P_ 

Substituting this into (19) and using 



(25) 



e i^S3P elsf 3Q = v^up) u^p) = exp [ ius p^ jspQ e ^/?P (26) 
which follows from (4), we obtain 



W(9) g^ s = exp [i(u 2 - s 2 )a(3/2 + iivs[3 2 ] e < auJ+ ^ Q e~^ s -^ p g\ (27) 

which is (20). Eq. (23) follows from Mehler's formula for the kernel of W(9) 
(cf. Merzbacher [5], p. 159, where 9 = —t). That g^l 2 = e t7r ^g follows from 
(23) by taking the limit 9 — > n/2. The expression for g n7T then follows from 
W(nn) = (W(n/2)) 2n and g(t) = g(-t). I 
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Corollary 2. Let T = {g u ,s '■ (u, s) G T} be any frame of Gabor wavelets, 
labeled by a discrete set V C R 2 , with window g and frame bounds < 
A < B < oo. For every real 9, let F = {(oj(9), s(6)) : (oj,s) G V} be 
the rotated version of V and let g e = W{9) g. Then the set of vectors 
Te = {g ' g : (a;, s) G Tq} forms a Gabor frame with window g and identical 
frame bounds A, B. 

Proof: For any (ou, s) G R 2 , Theorem 1 gives 

gl {e)Ae) =l{uj,s,9)- 1 W{9)g^ s . (28) 
Hence the Gabor transform with respect to the vectors in Te is 

( 9 i(O)XO) » / > = ^ s ' Q ) < W ^,s J)= 7(", «, f ) ( , ^M) / ) 
= -y(u>,s,0){g u , a ,f e ) 

(29) 

by the unitarity of W(9), where / 8 = W(-9)f G L 2 (R). Thus 

E K^^/)l 2 = E l(^iw,-w»/>l 2 = E KsW>l 2 - (so) 

(u,s)er 9 (w,«)er (w,«)er 

Since JF is a frame, 

4I/T< E KaW'>l 2 <s|l/T- (si) 

(w,«)er 

By the unitarity of W(-9), \\f e \\ 2 = \\f\\ 2 . Hence (30) shows that T e is a 
frame as well, with the same frame bounds A and B. I 

For example, if T is the frame generated by a compactly supported 
window g as in (10), then Te is generated by the window g e . This gives 
a one-parameter family of frames Te-, where only jF n7r have compactly sup- 
ported windows (n G Z). However, the frames Te with 9 = (n + |)7r have 
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windows which are compactly supported in the frequency domain. The fam- 
ily {Tq} therefore interpolates the two familiar types of easily constructed 
Gabor frames having windows with compact supports in time and in fre- 
quency, respectively. We call the family Tq a deformation of the frame T . 
Note that since g 2l7 (t) = —g(t), the frame jF 27r is essentially equivalent to 
JF = JF, and we only get distinct frames for < 9 < 2n. 

After this note was completed, I learned that work along similar (but 
inequivalent) lines has been done by R. G. Baraniuk and D. L. Jones [6]. They 
consider deformations of Gabor frames by operators that shear or chirp the 
time-frequency plane. I thank Bruno Torresani for pointing this work out to 
me. 
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